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ABSTRACT 

The orbital motion of a neutron star about its optical companion presents a window through 
which to study the orbital parameters of that binary system. This has been used extensively 
in the Milky Way to calculate these parameters for several high-mass X-ray binaries. Using 
several years of RXTE PCA data, we derive the orbital parameters of four Be/X-ray binary 
systems in the SMC, increasing the number of systems with orbital solutions by a factor 
of three. We find one new orbital period, confirm a second and discuss the parameters with 
comparison to the Galactic systems. Despite the low metallicity in the SMC, these binary 
systems sit amongst the Galactic distribution of orbital periods and eccentricities, suggesting 
that metallicity may not play an important role in the evolution of high-mass X-ray binary 
systems. A plot of orbital period against eccentricity shows that the supergiant. Be and low 
eccentricity OB transient systems occupy separate regions of the parameter space; akin to the 
separated regions on the Corbet diagram. Using a Spearman's rank correlation test, we also 
find a possible correlation between the two parameters. The mass functions, inclinations and 
orbital semimajor axes are derived for the SMC systems based on the binary parameters and 
the spectral classification of the optical counterpart. As a by-product of our work, we present 
a catalogue of the orbital parameters for every high-mass X-ray binary in the Galaxy and 
Magellanic Clouds for which they are known. 
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X-ray outbursts in transient high-mass X-ray binaries (HMXB) oc- 
cur when a compact object accretes intermittently from material 
being expelled by a massive companion. This intermittent accre- 
tion generally occurs because of an eccentric orbit, though tran- 
sient outbursts fro m near-circular orbit systems have been seen (e.g. 
IPfahl et al.ll2003) . Classically, such transient outbursts are thought 
to fall into one of two types. Type I outbursts occur near perias- 
tron passage of the neutron star (NS), where there is more acc- 
retable material. Studying these outbursts allows orbital periods to 
be derived through timing analysis of their light curves. Type II 
outbursts are typically an order of magnitude brighter than Type I 
outbursts and last much longer. The cause of this being an enlarge- 
ment of the circumstellar material or radially driven w ind around 
the companion star (e.g. iNegueruela & Okazakillioooh . allowing 
accretion to occur at any orbital phase and at a higher rate. Both 
of these transient outbursts are thought to originate from the Be/X- 
ray binary (BeXRB) class of HMXB in which the optical star is an 



early type Be or late Oe main sequence star and the compact ob- 
ject is a NS (there are currentl y no c onfirmed Be/black-hole bina- 
ries; iBelczvnski & Ziolkowskij2009i) . The NS is usually in a wide, 
eccentric orbit allowing for long periods of little X-ray activity. A 
third type of transient X-ray outburst has recently been detected and 
is now accepted to originate from a new class of HMXB known as 
supergiant fast X-ray transients (SFXT). These sources are char- 
acterised by rapid X -ray outbursts on the order of minutes to hours 
dSguera et al.l2005h . large dynamic ranges of order lO^'-lO^ and are 
associated with OB supergiant companion stars dChatv et al.l2008h . 
The other major sub-class of HMXB are the supergiant X-ray bina- 
ries (SGXB) in which a NS or black hole orbit a giant or supergiant 
star. Here the compact object is in a tight, near circular orbit, allow- 
ing for continuous accretion from the stellar wind or via Roche- 
Lobe overflow. The final sub-class making up the known HMXB 
family is the small, but growing class of y-ray binaries, in which 
very high energy radiation is thought to be produced via interaction 
of th e pulsar wind w ith the stellar atmosphere of the secondary star 
(e.g. lHiIletalj2010l) . 
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Table 1. Outburst date (To), duration and peak X-ray luminosity of the four 
systems investigated in this study. 



Source 


To (MJD) 


Duration (days) 


PeakL, (10"ergs s"') 


SXP2.37 


51573 


84 


21.0 


SXP6.85 


55435 


59 


3.3 


SXP8.80 


55080 


39 


7.3 


SXP74.7 


55232 


42 


3.5 



The Small Magellanic Cloud (SMC) holds a disproportion- 
ately large number of HMXBs when compared to the Milky Way. 
Based on their relative masses, there are a factor of 50 mo r e than 
would be expected. Population synthesis models by iDravl ( l2006l) 
predict the low metallicity environment in the SMC could increase 
the number of HMXBs by a factor of 3. However, this is still 
not sufficient to explain the SMC population on its own. A re- 
cent episode of star formation, possibly due to an increase in the 
tidal force exerted on the SMC by a close approach w ith the Large 
Magellanic Cloud (LMC: lOardiner & Noguchi"l996^. is the cur- 
rent favoured scenario to explain the number of binary systems 
seen. The metallicity difference is of particular relevance to our 
investigation as it has substantial influence in the evolution of mas- 
sive stars. The line-driven stellar wind s in massive stars are weaker 
in low metallicity environments (e.g. K udritzki. Pauldrach Sc Puis! 
fT987h . resulting in lower mass and angular momentum losses from 
the binary systems in which they reside. This could result in com- 
pact object masses different to those seen in the Galaxy and could 
even affect the evolutionary path of the binary system. The lat- 
ter may manifest itself as differences in observable parameters of 
the population, such as spectral distribution, the distribution of or- 
bital sizes and eccent ricities or the distribution of mass functions. 
iMcBride et feOOSi) show that, despite the low metallicity in the 
SMC, the distribution of spectral types of the optical counterparts 
to Be/X-ray binaries (BeXRB) is consistent with that of Galactic 
systems. 

Our aim is to use several years of data taken with the Rossi X- 
ray Timing Explorer (RXTE) observatory to uncover the first sizable 
sample of orbital parameters of HMXBs in the SMC. We do this by 
fitting a radial velocity model to NS spin period variations that are 
seen during several Type 11 outbursts. If the variability is due to the 
motion of the NS around the optical counterpart, we can obtain a 
best fit solution of the orbital parameters. The model and data are 
presented in section 2. The nature of this method means a good 
fit can only be achieved if the outburst lasts for at least one orbit 
and the observations are made regularly throughout that orbit. In 
section 3 we sho w the fits to 4 syste ms that were possible with our 
data sample (see lCalache et al.' 2008 for details of our monitoring 
programme). These results now make 6 systems in the SMC for 
which we have good orbital solutions (excluding the only known 
SGXB SMC X-1). In section 4 we discuss these systems in the 
context of the Galactic population before summarising our results 
in section 5. 



2 DATA & ORBITAL MODEL 

RXTE monitoring of the SMC dCalache et ai] l2008^ has resulted 
in 10 year X-ray light curves and spin period histories for most of 
the known pulsars in the SMC. During this time, we have detected 
12 systems that have shown evidence for one or more Type II out- 
bursts. Of these systems only a handful were of sufficient length and 
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Figure 1. The spin period of SXP6.85 as measured by RXTE during 2 
months of pointed observations. Both the orbital motion of the NS and the 
accretion driven spin-up are visible. Overplotted is the model fit to the data 
as shown in Table |2] 



had frequent enough observations within the outburst to achieve a 
good model fit. Table[T]gives a summary of the start time, duration 
and peak X-ray luminosity in the 3-10 keV band, of the studied out- 
burst for each system. Each observation was made using the RXTE 
PCA in GoodXenon mode. Data were extracted in the 3-10 keV 
energy range and binned at 0.01s. The resulting light curves were 
background subtracted and barycentre corrected before finally cor- 
recting the count rate for number of active PCU's. The final light 
curves were passed through a Lomb-Scargle periodogram to search 
for pulsations. The error associated with any period found is cal- 
culated based on the formula for the standard deviation of the fre- 
quency given in .Home & Baliunas, (il986i) : 



(1) 



where cr^ ^ is the variance of the light curve, N is the number of 
data points, D is the length of the data and A is the amplitude of the 
signal given by: 



N 



(2) 



where zo is the Lomb-Scargle power. 

In the simplest example of the orbital plane being perpendic- 
ular to our line-of-sight (inclination, i = 0), any change seen in the 
spin period of the NS should be due to torque produced by accretion 
onto the NS surface. The most common result of this is for the spin 
period to decrease over the length of the outburst. This is known 
as spin-up. However, in some cases it is seen that the spin period 
increases during phases of accretion. This is known as spin-down. 
Spin-down also occurs during phases in which there is no accre- 
tion due to torques exerted by the magnetic field of the NS. In the 
more generic case of the orbital plane being at some angle to our 
line-of-sight, we see more complicated spin period variations. Con- 
volved with the standard spin-up of the NS, there will be a shifting 
of the X-ray pulse arrival time caused by the orbital motion of the 
NS around its counterpart. It is because of this that our orbit fit- 
ting code simultaneously fits a simple spin-up model and an orbital 
radial velocity model 
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Table 2. The orbital parameters for SXP6.85 from the analysis of 3-10 keV 
RXTE PCA data. 



Parameter 



Orbital Solution 



Orbital period 

Projected semimajor axis 

Longitude of periastron 

Eccentricity 

Orbital epoch 

Spin period 

First derivative of P 

Goodness of fit 



Porbital (d) 

a^&mi (light-s) 
e 

'^periastron (MJD) 

P(s) 
PClO-'Oss-') 

Xv 



21.9 ±0.1 
151 ±6 
125 ±6 
0.26 ± 0.03 
55479.1 ±0.4 
6.8508 ± 0.0001 
-8.0 ±0.5 
1.79 



1 + 



v(0 



P{t) 



(3) 



where P„i,s is the detected spin period of the pulsar, c is the speed 
of light and v{t) is the binary radial velocity. v(f) is calculated us- 
ing the IDL routine binradve l.pro from the aitli^B which uses the 
proced ure of iHilditcl] ( 1200 ih and employs the method of iMikkolal 
( Il987h to solve Kepler's equation. The spin of the pulsar is given 
by P{t) which incorporates the spin-up component and is defined 
below: 



P(t) 



P{to) + Pit - to) - 



P{t-taf 



(4) 



where P and P are the spin-up and change in spin-up of the NS. 
We perform a least-squares fit on the data and return the best fit- 
ting parameters. The full spin-up component shown in equation |4] 
was fit to every outburst. However, no data set allowed for suffi- 
cient statistics in the fit to justify a P component and hence it was 
removed. 

The results of our fits are presented in the next section. In order 
to obtain a meaningful fit, the period measurements used in the fit- 
ting procedure were restricted to detections greater than a particular 
significance. This threshold was not kept constant for every source 
because the luminosity and duration of each outburst was different, 
meaning fits improved or worsened by including data above dif- 
ferent thresholds. As such, a significance threshold was chosen to 
minimise the (Icr) errors in the binary parameters for each fit. The 
chosen value is stated in the relevant section. 



3 RESULTS 

3.1 SXP6.85 = XTE J0103-728 

SXP6.85 was first detected in 20 03 by RXTE as a 6.848 s pulsed 
X-ray source dCorbet et al.ll2003l) . It was later detected in a 2006 
XMM-Newton observation at the position R.A. = 0l''02'"53U, dec. 
= -72°44'33'.'0 (J2000.0). This d etection led to the identification 
of a V=I4.6 optical counterpart (lHaberletal.ll2008h, allowing it 
to be classified as an HMXB. Follow-up work by l lVlcBride et al.l 
( l2008h classified the counterpart as an 09.5 V-BO V emission line 
star. In subsequent years it has been detected on 5 distinct occa- 
sions, coinciding with times when the counterpart was optically 



' Developed at the Institut fiir Astronomic und Astrophysik, Abteilung As- 
tronomic, of the University of Tubingen and at the University of Warwick, 
UK; see http://astro.uni-tuebingen.de/software/idl/ 
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Figure 2. The spin period of SXP2.37 as measured by RXTE. The orbital 
modulation is not as obvious as in the case of SXP6.85, but the period de- 
viations around the general spin-up trend suggests we are seeing orbital 
modulations. Overplotted is the model fit to the data as shown in Table [3] 



bright jTownsend etal]|20I 0'). 'McGo wan eTal] ([20081) find that the 
optical flux varies by ~ 0.5 magnitudes with a period of 620 ± 18 d. 
They associate this with the growth and decay of the circumstellar 
disk. Those authors also show that the source gets redder as it gets 
brighter, suggestive of a low inclination system, and propose a low 
eccentricity based on comparison to other systems. 

Until now, the orbital period of this sy stem was not known for 
cert ain. Analysis of optica l light curves bv lMcGowan et al" I( l2008l) 
and lSchmidtke & Cowlevi ( |2007|) showed hints of periodicities at 
1 14. 1 ± 0.6 d and 24.8 ± 0. 1 d, but neither could be confirmed as the 
orbital period of the system. The 25 d period was c loser to the ex - 
pected orbital period based on the Corbet diagram (ICorbet,lll98^ . 
but the 1 14 d period was reinf orced by the detectio n of a 1 12 d pe- 
riod in the X-ray light curve dGalache et al.ll2008! >. To try and re- 
solve this issue, the two longest Type II X-ray outbursts were fit 
with our orbital model. The outburst occuring around MJD 54800 
turned out to be too sparsely covered, meaning we were unable 
to get an acceptable fit. However, the outburst beginning on MJD 
55435 (Table[T]l was sampled very well thanks to dedicated point- 
ings at the position of the source in addition to our regular monitor- 
ing. The period evolution is shown in Fig. [T] with the best model 
fit overplotted. The data used in the fit were cut at the 99.99 % 
(4cr) significance level. The very clear changes in the spin period 
allowed the radial velocity of the NS to be found quite simply. Un- 
fortunately, we were unable to fit both outbursts simultaneously as 
the amount of spin-down happening between the outbursts is un- 
known. The solution is presented in Table |2] and is the best fit we 
have in our sample. We propose 21.9 ± 0.1 d to be the true orbital 
period of this system, placing it nicely in the BeXRB region of the 
Corbet diagram. The other orbital parameters and the nature of the 
25 d and 1 14 d periods are discussed in section 4. 



3.2 SXP2.37 = SMC X-2 

SXP2.37 was d iscovered by SAS 3 observations in 1977 
I Li. Jemigan & Cla rk 1977) as a highly variable X-ray source. 
ICorbet et al.l ( 1200 ih discovered 2.372 s pulsations from the source 
during the Type II outburst from January 2000 to May 2000. The 
proposed optical counterpart to this X-ray source was shown to be 
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Table 3. The orbital parameters for SXP2.37 from the analysis of 3-10 keV 
RXTE PCA data. 



Parameter 



Orbital Solution 



Orbital period 

Projected semimajor axis 

Longitude of periastron 

Eccentricity 

Orbital epoch 

Spin period 

First derivative of P 

Goodness of fit 



Pm-bital (d) 

OxSmi (light-s) 
e 

'^periaslron (MJD) 
P(S) 

/'(10-"'ss-l) 



18.38 ±0.02 
73.7 ±0.9 
226 ±8 
0.07 ± 0.02 
51616.8 ±0.4 
2.37194 ±0.00001 
-0.720 ± 0.015 
3.58 



a close north-south double ( iMurdin. Morton & Thom"aslll979l) . As 
such, it has been difficult to tell which of the two stars is the true 
optical counterpart due to large X-ray error circles and the lack 
of any significa nt periodicity in OGLE II and the fir st 5 years of 
OGLE III data jSchmidtke, Cowley & UdZskil 12006 ). It was not 
until a further 2 years of OGLE HI data became availa ble that the 
correct counterpart revealed itself . Ischurch et allteOl ll) show there 
is a clear 18.62 ± 0.02 d periodicity in the final 2 years of the light 
curve of the northern star, making it very probable that this star is 
the true counterpart and that the periodicity is the orbital period of 
the binary system. ISchmidtke. Cowley & Udalskil ( |2009|) suggest 
that the 18.6 d period, along with the apparent 9 and 6d harmon- 
ics a re caused by 2 nonradial pulsations seen around 0.86 s and 
0.9s. lSchurch et al.l (201 1) argue that because the 18.6 d modula- 
tion is visible in the light curve and the folded light curve at this 
period is typical of other binary system s, the modulation must be 
orbital in nature. iMcBride et al.l ( l2008h classify the counterpart as 
an 09.5 III-V emission line star. 

The outburst that led to the discovery of pulsations from 
SXP2.37 ( ICorbetetal1l200lh is the only time RXTE has detected 
the source with any high significance or for any great period of 
time. Thus, we tried fitting the period measurements from the out- 
burst with our model. Fits were tried both with the orbital perio d 
free to vary and fixed at the value reported in ISchurch et allfcoilh . 
Both fits yielded the same parameters to within their errors and very 
similar goodness of fits, verifying that the 18.6d periodicity is the 
orbital period of the system. Fig.|2]shows the pulse period evolution 
of the source. The data are cut at the 99.9999% (5 cr) significance 
level. Overplotted is the best model fit to the data in which the or- 
bital period was free to vary. We refine the orbital period measure- 
ment to 18.38 ± 0.02 d. The other parameters are shown in Table|3] 
and are discussed later. 



3.3 SXP8.80 = RX J0051.8-7231 = 2E 0050.1-7247 

SXP8.80 was discover ed by Einstein IPC observations in 1980 
teruhweiler et al.ll 19871) as a new X-ray source in the SMC wing. 
The original error circle from these observations l ed to a ten- 
tative association with the optical star AzV 111. 'Wang & Wu' 
( Il992l) showed that this source was highly variable, had a mod- 
erate X-ray flux and a hard spectrum, suggesting a BeXRB clas- 
sification. ROSAT PSPC observations in May 1993 discovered a 
periodicity of 8.8 8163 ± 0.00001 s in the direction of SXP8.80 
dlsrael et al. Those authors also showed the transient na- 

ture of the source which, along with the spin period, confirmed 
this as a new BeXRB in the SMC. A refined X-ray position pub- 
lished bv lKahabka & PietschI (199^ cast doubt over the true opti- 
cal counterpart, with more than one star looking possible. It was 
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Figure 3. The spin period of SXP8.80 as measured by RXTE. The orbital 
modulation is apparent, though the data are sparse over some of the orbit. 
The orbital period had to be frozen at 28.51 d to produce an acceptable fit. 
This fit is overplotted and is shown in Table|4] 



not until i Haberi & Sas^ 1200(f) improved the ROSAT error cir- 
cle that the star MA93 506 was identified as the most probable 
counter part. The orbital pe riod of SXP8.80 was found to be 28.0 
+ 0.3 d ('Corbet et al.ll2004l) during RXTE monitoring of the SMC. 
Subsequent analysis of MACHO and OGLE data of MA93 506 re- 
vealed periodicities of 185 d and 33 d respectively dCoe et alj2005l: 
.Schmidtke & Cowley 2006), casting doubt over the orbital period 
and the optical counterpart determination. iMcBride et al.l ( 1200 8h 
classify MA93 506 as a 09.5-BOIV-V emission line star. Only 
recently has this star been verified as the counterpart with the de- 
tection of a 28. 51 ± O.OI d periodicity in the com bined OGLE II & 
III light curve dRaioelimanana. Charles & Udalskii,201 1,) . in agree- 
ment with the X-ray derived period. 

SXP8.80 had been X-ray quiet since a long period of Type 
I outbursts that began around MJD 52700. The Type II outburst 
to which we have fit our model began on MJD 55080 and is one 
of the most luminous outbursts seen in the SMC (see Table [T); 6 
or more harmonics of the NS spin period were often seen in the 
power spectra. Unfortunately, the coverage during this outburst was 
sparse and so there are large gaps throughout the outburst. The data 
are shown in Fig. |3] and are cut at the 90 % significance level to 
try and include as much data as possible. This means that some 
of the data toward the end of the outburst have large associated 
errors. To ai d fitting further, the orbital period was frozen to 28.51 d 
dRaioelinianana. Charles & Udalski,2011i) . Our best fit is plotted 
over the data in Fig. |3] The results are given in Table |4] and are 
discussed in section 4. 



3.4 SXP74.7 = RX J0049.1-7250 = AX J0049-729 

RX J0049. 1-7250 was discovered d uring a ROSAT PSPC ob ser- 
vation of the SMC in October 1991 ( Kahabka & Pietschlll99d) . It 
appeare d variable by more t han a factor of 10 and was highly ab- 
sorbed. iKahabka & PietschI jl99d) concluded that the source was 
an XRB on the far side of the SMC, but could not rule out a 
time variable AGN. A new pulsar was discovered by RXTE dur- 
ing p ointed observations of SMC X-3, with a periodicity of 74.8 ± 
0.4 s dCorbet et al.|[T998h . This was confirmed by the detection of 
a 74.675 ± 0.006 s period during an ASCA observation on 1997 
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Table 4. The orbital parameters for SXP8.80 from the analysis of 3-10 keV 
RXTE PCA data. 



Parameter 



Orbital Solution 



Orbital period 

Projected semimajor axis 

Longitude of periastron 

Eccentricity 

Orbital epoch 

Spin period 

First derivative of P 

Goodness of fit 



Porbital (d) 

flvsinf (light-s) 
e 

'^periaslroii (MJD) 
P(S) 

PClO-'Oss-') 

1 

Xv 



28.51 (frozen)^ 
112 ±5 
178 ±4 
0.41 ± 0.04 
55106.7 ±0.2 
8.9038 ± 0.0001 
-6.9 ± 0.6 
2.60 



* Orbital period frozen at the period found by 
iRaioelimanana. Charles &UdalskiH201lh 



November 13 llYokogawa & Kovamj [1991) . iKahabka & PietschI 
( Il998h associated th e A5CA source with the highly variable ROSAT 
source JiCahabka & PietschI [l996) and suggested this as a Be 
type tr ansient. Optical follow up work bv lStevens. Coe & Buckley! 
( Il999h deduced one probable and one possible counterpart to 
the X-ray source (which they refer to as objects 1 & 2). Ob- 
ject 1 has since been confirmed as the correct counterpart from 
the detection of a 33.4 ± .4 d period in the optical light curve 
dSchmidtke & CowlevlllOOSi). assumed to be the orbital period of 
the svstem. lMcBride et alJiioOSi) classify the counterpart as a B3 V 
emission line star. 

There have been sporadic detections of SXP74.7 with RXTE 
since its discovery, although no major outburst was seen until MJD 
55232. This Type II outburst was observed regularly by RXTE, to 
try and follow the period evolution. Unfortunately, the outburst 
lasted for less than one orbit, making fitting difficult. The period 
history is plotted in Fig.E) with the data cut at the 99.9999% (5 tr) 
significance level. Besides a small kink at around MJD 55247, there 
is very little binary motion visible making fitting more challenging. 
This could be caused by a low orbital inclination with respect to our 
line-of-sight or the modulation could be getting swamped by an ex- 
ceptionally large spin -up component. The orbital period was froz en 
to the value found bv lRaioelimanana. Charles & Udalskil | |201 ih to 
help fit the data. Our best fit is plotted over the data and the pa- 
rameters are given in Table |5] From the reduced chi-square value 
it is clear that the model is over-fitting the data. This is reflected in 
the larger errors on many of the parameters. Again, the results are 
discussed in the next section. 
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Figure 4. The spin period of SXP74.7 as measured by RXTE. Very little 
orbital modulation is apparent, though the outburst is sampled well. To aid 
fitting, the orbital period was frozen to 33.38 d. The fit is overplotted and is 
shown in Table [5] 



Table 5. The orbital parameters for SXP74.7 from the analysis of 3-10 keV 
RXTE PCA data. 



Parameter 



Orbital Solution 



Orbital period 
Projected semimajor axis 
Longitude of periastron 
Eccentricity 
Orbital epoch 

Spin period (MJD 55247.2) 
First derivative of P 
Goodness of fit 



Porbilal (d) 

flvsini (light-s) 
a>n 
e 

'^periastron (MJD) 

P(s) 
/■(lO^'^ss-') 



33.38 f/rozenj' 
147 ± 15 
186 ± 18 
0.40 ± 0.23 
55280.9 ± 1.5 
74.867 ± 0.002 
-635 ± 37 
0.33 



' Orbital period frozen at the period found by 
IRaioelimanana, Charles & Udalski 1 2011.) 



.fx(M) -. 
and 



K/P{\-e^) 
2nG 



Mc sin^ i 



(5) 



4 DISCUSSION 

The 4 systems presented above bring the total number of BeXRBs 
with reliably measured orbital parameters in the SMC to 6. As such, 
it is the first time in which we have a large enough sample of SMC 
systems to compare with the parameters of Galactic HMXBs. In 
this section we derive mass functions, inclination angles and or- 
bital semimajor axes for our SMC sample and compare them and 
the other binary parameters to parameters calculated from studying 
Galactic systems. 



4.1 The Binary Mass Functions 

The masses of the two stars can be described by the X-ray and 
optical mass functions of the binary system: 



/c(M) 



/<:,V(1 -^2)3/2 M,.sin'/ 



2;rG 



(1 + 1/?)- 



(6) 



respectively. and are the masses of the NS and the optical 
counterpart, / is the inclination of the orbital plane to the line of 
sight, P is the orbital period, e is the eccentricity and q is the mass 
ratio (=M^/Mc). The semiamplitudes of the radial velocity curves 
are given by 



K, 



2na„ sin ; 



(7) 



/'(I -e2)i/2 

where a,, and a,, are the semimajor axes of the ellipse travelled by 
the NS and the optical counterpart about the centre of mass of the 
system. Thus, if the radial velocity curves for both stars are known, 
along with the inclination angle, one can calculate the masses of 
the two stars to high precision. However, due to the non-eclipsing 
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Table 6. The derived X-ray mass function, orbital inclination and semi- 
major axis for the 4 systems studied in this paper and the 2 previously 
studied SMC BeXRBs. is estimated from the spectral classification 
of the authors referenced here using the table of luminosity and spec- 
tral type from de Jager & Nieuwenhuiizen ( 1987), along with the standard 
mass-luminosity relationship for main sequence stars (L oc M^ *). fx(M) is 
calculated using results from our fits and is shown in solar masses. The 
inclination and semimajor axis are estimated using equation |5] and 
M, = 1.4Mo. 



Source 


(M,/Mo) 




in 


a, (R*) 


SXP2.37' 


18-26 


1.27 ±0.05 


24 +\ 


9 


SXP6.85' 


13-26 


7.71 ±0.92 


50 !f 


12 


SXP8.80' 


13-23 


1.86 ±0.25 


30 


13 


SXP74.7' 


6-9 


3.06 ± 0.94 


55 +^0 

-13 


16 


SXPll.5^ 


13-26 


3.79 ± 0.48 


37!" 


16 


SXPI8.33 


8-13 


1.43 ±0,17 


36 « 


10 



Note. R eferences for the spectra l classification a r e give n: iMcBride et al.l 
i2008l) , frownsend et al." ( 201 1), ^ Schurch etaP j2009 l). The spectral type 
of SXP18.3 in Schurch et al. (2009) is estimated based on optical an NIR 
photometry as no spectrum allowing for classification is available, 

nature of these SMC binary systems, we are unable to make pre- 
cise measurements of the NS masses. Instead, we can calculate the 
X-ray mass function for each of our systems using our orbital solu- 
tions and make estimates of the inclination and orbital size a, ) 
using masses estimated from the spectral classification of the coun- 
terpart and the standard mass of a NS of 1.4 Mq. These results are 
shown in Table |6] for the 6 SMC BeXRBs with orbital solutions. 
The range in the value of each inclination is based mostly on an 
uncertainty in the spectral classification of one spectral type either 
side of the published value. The uncertainty in the mass function 
also contributes to the range in inclination, albeit less significantly. 
The semimajor axes are estimated based on the most probable mass 
and radius of the primary star. 

The orbital inclinations seen in this sample are as expected 
given the method of detection. Very low inclinations would mean 
the delays in pulse arrival times would not be detected, whilst very 
high inclinations would mean the X-ray source gets eclipsed by the 
primary star. More specifically, we can compare the inclinations 
to Ha profiles of the Be star to investigate the (mis-)alignment of 
the orbital plane and the circumstellar disk. We find that SXP6,85 
and SXP74,7 have quite narrow, single peaked Ha emission, de- 
spite having the highest estimated inclinations, Conversly, we note 
that SXP18,3 has a lower inclination but shows prominent dou- 
ble peaked Ha emission in its optical spectrum. However, we note 
that the mass of SXP18,3 has been estimated photometrically, not 
spectrally. Although this may provide some qualitative evidence 
of orbit-disk misalignment in these systems, more quantitative ev- 
idence may only be obtained from detailed polarimetric studies 
of the disk itself. The estimated semimajor axes in Table |6] are 
quite similar and range from 9 to 16 stellar radii, but are these 
values what one might expect? The use of complex models that 
describe the dynamics and thermal structure of the circumstellar 
disk around the primary star helps explain what is observed here, 
lOkazakil j2007h shows that the density of the disk is several orders 
of magnitude lo wer beyond 10 R* than it is closer to the star, whilst 
ICarciofil ( 1201 Ol) shows most of the optical and NIR flux and po- 
larisation is emitted from within 10 R* and nearly all FIR and Ha 
flux is produced inside 20 R*, Thus, these models are predicting 



that a very large fraction of the matter in the disk is within ap- 
proximately 10 R*, This prediction goes some way to explaining 
the nature of the X-ray outbursts that we see in the SMC binary 
systems, SXP74,7 and SXP11,5 have the largest predicted orbital 
size (relative to the radius of the counterpart) and are observation- 
ally the least active of the sample, rarely undergoing a Type I out- 
burst and only once being seen in a Type II outburst, SXP6,85, 
SXP8,80 and SXP18.3 have smaller predicted orbital sizes and are 
much more active, often being detected in Type I or Type II out- 
bursts. This is likely due to the NS passing closer to, or further 
into, the circumstellar material than those in larger orbits. The ex- 
ception to this seems to be SXP2,37 which has the smallest pre- 
dicted orbital size, but is very rarely detected in X-ray outburst. 
This can be explained when conside ring the low eccentricity of the 
orbit, lOkazaki & Negueruelal (1200 ih predict that for low eccentric- 
ity HMXBs, the circumstellar disk will get truncated at the 3:1 res- 
onance radius, in contrast to intermediate and high eccentricity sys- 
tems in which truncation is much less efficient and occurs at larger 
distances near the Roche lobe radius of the star Thus, the disk in 
SXP2,37 could be truncated at a much smaller radius than those of 
the other systems in our sample, explaining the small number of 
X-ray outbursts seen, 

4.2 Comparision with the Galactic population 

In order to properly compare this new sample to those found in 
the Galaxy, a thorough review of the literature was made and a 
database compiled of all known HMXB systems for which the or- 
bital period and the eccentricity are known. This list is intended 
to be exhaustive, although this cannot be guaranteed. In total we 
find 25 HMXB systems with a Be companion and 17 with a gi- 
ant or supergiant companion in addition to the 6 systems discussed 
so far in this paper. Table |7]presents all 48 HMXBs with their or- 
bital period and eccentricity, as well as spectral type, spin period, 
projected semimajor axis, longitude of periastron and mass func- 
tion where known. O ur source list i s based mostly, but not so lely, 
on tables presented inlsildsten et all ( ll997h , |Pfahl et alj ( |2002|) and 
iMartin, Tout & P ringle ('2009^, We also used the updated database 
of BeXRBs of [Raguzova & Popov (2005) and the online database 
of INTEGRAL sources maintained by Jerome Rodriguez and Arash 
Bodaghe^B In Fig, [5] we plot the orbital and spin periods of these 
HMXBs (less 9 systems with unknown spin periods or PSR des- 
ignations). The red diamonds represent the SG systems, the blue 
triangles are the Galactic BeXRBs and the green stars are the SMC 
BeXRBs, As ex pected, the SG and Be systems occupy separate 
parts of the plot Jcorbet.ll 19861) and the SMC BeXRBs sit amongst 
the Galactic Be systems. This sample will be used to compare with 
the results from our analysis throughout the discussion. 

An observation that emerges from this work, is that there are 
often many different periodicities found in the optical and X-ray 
light curves of individual sources and that these periods do not 
always agree. As such, great care must be taken when deciding 
which periodicity is the orbital period and which can be attributed 
to stellar pulsations, disk oscillations, superorbital periods or even 
aliasing between different periodicities. Fig, |5] supports our inter- 
pretation of the 21,9 d orbital period in SXP6.85 as a longer 1 14 d 
period would see it sit towar ds the edge of th e BeXRB region. 
The 1 12 d period reported bv lGalache et al,l ( |2008 ) came from tim- 
ing analysis of the long term RXTE light curve, A periodogram of 

^ http://irfu.cea.fr/Sap/lGR-Sources/ 
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Table 7. All high-mass X-ray binaries with reliably known orbital periods and eccentricities. 



Source name 


Spectral type 


P spin (s) 


Porhital (d) 


e 


c/vsin/ (light-s) 




fxiM) (Mo) 


Be/X-ray binaries 
















r5>K JUU4j— /3iy 


Bl V 


0.93 


c 1 1 An 
Jl.loy 


n ono 


\n A TIC _i_ A A7^T 

1 /4.Z3j + U.UOZ 


1 IC TIA _i_A A AT 
1 1 J.Z36 + O.OUZ 


2.169 


KA JUU4y. i— iIjU (oAr /4. / j 


153 V 


1 A Q 

/4.y 


11 iQ _i. n ni 
33. 3o ± U.Ul 


A APi _i. A 01 

U.4U + U.Z3 


1 yl"? _i. 1 C 
14/ + 15 


1 OA j_ 1 C 

I06 + lo 


1 nA -1. A Q/1 

3.U0 + U.y4 


Zr, UUjU. 1— /24/ (oAro.oUJ 


wy.j-xiUlV-V 


o.oy 


CI _i. n ni 
Zo. J 1 ± U.Ul 


n /1 1 -1. n f\A 
U.41 + U.U4 


1 1 T _i. C 
1 IZ + J 


1 TO j_ A 
1 /O + 4 


1 OA -1. A TC 
1 .00 + U.Zj 


MA UUj J+OU4 (y L.as ) 


BO. 5 IV 




oni cn _i_ n Tn 

zo3.3y ± o.zy 


n lAn _i_ n mc 
U.ZoU + U.U3j 




An n _t_ Q n 
4/.y + 0.0 




Air, JuUjj— /Z / (iSArio.jJ 


B 1-B3 V 


18.4 


1 1 ^n _i_ n n 1 
1 /. /y ± U.Ul 


n /1 1 _i_ n m 
U.43 + U.U3 


75 + 3 


1 C _i_ A 

1 J + 6 


1 .43 + 0. 17 


iMC A— Z (oArZ. J / J 


(jy.j 111- V 


2.37 


1 Q 10 _i_ n m 
lo.3o ± U.UZ 


n m _i_ n m 
U.U/ + U.UZ 


11 "7 _i_ n A 

/3. / + o.y 


226 + 8 


1 T7 _i_ A AC 

1 .Z / + U.Uj 


AlnJUlUj— /Zo (oAro.ojJ 


(jy.j-iiU IV-V 


A oc 
O.SJ 


21.9 ±0.1 


n TA _i_ n m 
U.Zo + U.U3 


151 + 6 


125+6 


"7 "7 1 _i_ A AT 

/. / 1 + u.yz 


ICjK JU1054— 7253 (aXrllo) 


uy.j-Boiv-v 


1 1.5 


36.3 ± 0.4 


U.Zo + 0.03 


167 + 7 


224 + 10 


1 1(\ 1 f\ AO 

3. /y + 0.4o 


/IT I r\ 1 1 c 1 £;i5 

4U 0115+63 


BO. 2 V 


3.61 


24.317 


A 1 /1 1 1 r\ A A A 

0.34Z + 0.004 


1 /I A Ar\ 1 A "71 

140.69 + 0.7Z 


/IOC 1 f\ Cl 

4o.5 + O.y 


5.06 + O.Oo 


Lb 1 +ol 3U3 


BO V 




Zo.4y6 ± 0.003 


U.53 / + 0.034 




40.5 + 5.7 




VU33Z+J3 (l3ljL.am) 


08.5 V 


4.38 


lA cn _i_ n on 

3O.50 ± u.zy 


A /I 1 "7 _i_ A Am 

U.41 / + U.UO/ 


OT /I A _i_ A C\A 

oz.4y + u.y4 


TCI AC\ _i_ n ni 
Zo3.4y + O.y 1 


n /I c _i_ n nT 
0.4j + O.UZ 


4U U3jZ+3Z (AFer) 


(~\C\ Z. TTT TiCWT 

L>y.5 III-dO V 


837.7 


ocn 1 _i_ n A 
ZD0.3 ± U.O 


n 1 1 1 _i_ n ni o 
U. 1 1 1 + U.Ulo 


454 + 4 


TOO _i_ n 

Zoo + y 


1 An _i_ n r\A 
1 .00 + 0.04 


CI Cam^ 


B4 III- V 




iy.41 ± U.UZ 


n AT _i_ n n 1 
U.oZ + U.Ul 


140 + 4 




"7 OT _i_ n AT 

/.oZ + 0.0 / 


A (j53D+Zo 


BO III-V 


103.5 


1 1 n Q _i_ n "3 
1 10.3 ± 0.3 


n /I "7 _i_ n m 
U.4 / + U.UZ 


267 + 13 


130 + 5 


1 .68 + 0.25 


A (J53o— OD 


B2 III 


0.07 


1 A A/1 A 


n OT _i_ n (\A 
U.oZ + 0.04 




TTT _i_ T 1 
ZZZ ± Zl 




oAA JUo35.Z+Uj33 


TJ 1 TTT Ti'^ \T 

ol lll-iiZ V 


U.U3 


11 1 _i. n c 
1 1.Z ± O.J 


A TO _i. A AO 

O.zy + u.uy 


01 j_ 1 1 
o3 + 1 1 


1CA j_ T/1 

356 + Z4 


/I j_ T A 

4.y + z.o 


Ljo Uo34— 43U 


Tt(\ 1 TTT \T 

oU-Z 111- V 


1Z.3 


1 nc _i_ n /I 
IUj.o ± U.4 


A 1 /I _i. A A/1 

U. 14 + U.U4 


1 TO -1. /1 1 

IZo + 43 


1 /I A j_ A A 

140 + 44 


A T j_ A T 

u.z + u.z 


DCl? 1TCQ Aq14 
rjK IZjy— 03 


T30 


n nc 


1 TIA "70/1 

iZ30. /Z4 


A 0*7 

U.O / 




110 AAC 




(jKw J iUUo— J / 


B1-B2 V 


93.5 


0/1 "7 _i_ n /I 
Z4 /.o ± U.4 


A AO _i_ A AT 
U.OO + U.UZ 


530 + 60 




2.6 + 0.9 


1A 1110 /;i/;16 
lA 1 1 lo— OlO 


f~\C\ C T\7 A 7 

Kjy.j IV- V 


407.7 


T/1 n _i_ n /I 
Z4.U ± 0.4 


< 0.16 


54.9 + 1.4 


310 + 30 


All _i_ A Al 

0.3 1 + 0.U3 


TC 1 /I 1 "7 AT/1 17 

Zo 141 /— oZ4 


Bl V 


17.5 


/ii 1 T _i_ n m 
4Z. IZ ± U.U3 


A /I 1 "7 _i_ A AA1 
U.41 / + U.U03 


Tm 1 _i_ 1 n 
ZU/.l + l.U 


Tno oc _i_ n AO 
Zyo.OJ + U.OO 


c n _i_ A AO 

J. 3 / + U.Uo 


XTEJ1543-568^^ 


Be 


27.1 


75.56 + 0.25 


< 0.03 


353 + 8 




8.27 + 0.56 


2S 1553-542^'^ 


Be 


9.27 


30.6 ± 2.2 


<0.09 


164 + 22 




5.1 +2.2 


SWIFT J1626.6-5 1562** 


Be 


15.4 


132.89 + 0.03 


0.08 + 0.01 


401 + 5 


340 + 9 


3.92 + 0.15 


GROJ1750-272i 




4.45 


29.806 + 0.001 


0.360 + 0.002 


101.8 + 0.5 


206.3 + 0.3 


1.27 + 0.02 


2S 1845-024^2 


Be 


94.3 


242.18 + 0.01 


0.879 + 0.005 


689 + 38 


252.2 ± 9.4 


6.0+ 1.0 


4U 1901+0323 


OB 


2.76 


22.583 


0.036 


106.989 + 0.015 


268.812 + 0.003 


2.578 + 0.001 


GRO J 1944+2624 


BO-1 IV- V 


15.8 


169.2 + 0.9 


0.33 + 0.05 


640 + 120 


269 + 23 


9.8 + 5.5 


GROJ1948+3225 


BOV 


18.7 


40.415 + 0.010 


0.033 + 0.013 


137 + 3 


33 + 3 


1.69 + 0.11 


EXO 2030+37526 


BOV 


42 


46.021 


0.412 + 0.001 


244 + 2 


211.3 + 0.3 


7.36 + 0.18 


SAX 12103.5+45452"^ 


BOV 


358.6 


12.665 


0.406 + 0.004 


74.07 + 0.86 


244.3 + 6.0 


2.72 + 0.09 



Giant & Supergiant X-ray binaries 



SMCX-128 


BO lb 


0.72 


3.892 


0.0002 


53.577 


317±9 


10.897 


28 0114+6502' 


Blla 




11.598 


0.18 ±0.05 




51 ± 17 




LMCX-4^'' 


08 III 


13.5 


1.408 


0.006 ± 0.002 


26.343 ± 0.002 




9.893 ± 0.002 


VelaX-l^' 


B0.51a 


283.5 


8.964 


0.090 ± 0.001 


113.89 ±0.13 


152.59 ± 0.92 


19.73 ± 0.07 


CenX-332 


06-8 III 


4.8 


2.087 


< 0.0001 


39.661 




15.373 ± 0.001 


IE 1145.1-6141'^ 


B21a 


297 


14.365 ± 0.002 


0.20 ± 0.03 


99.4 ± 1.8 


308 ±8 


5.11 ±0.28 


0X301-23* 


Blla 


685 


41.498 ±0.002 


0.462 ± 0.014 


368.3 ± 3.7 


310.4 ± 1.4 


31.14 ±0.94 


4U 1538-5235 


BO lab 


526.8 


3.723 


0.174 ±0.015 


56.6 ±0.7 


64 ±9 


14.04 ± 0.52 


IGRJ16393-46433* 




910.4 


3.688 


0.15 ±0.05 


55 ±2 




13.1 ± 1.5 


IGR J 16493-4348" 


B0.51b 




6.782 ± 0.002 


<0.15 


< 128 




<49 


OAO 1657-41538 


B3Ia-ab 


38.2 


10.444 ± 0.004 


0.104 ±0.005 


106.0 ± 0.5 


93 ±5 


11.72 ±0.17 


EXO 1722-3633' 


BO-1 la 


414.8 


9.740 


<0.19 


101 ±3 




11.7 ± 1.0 


PSR J 1740-3052'"' 


B 


0.57 


231.030 


0.579 


756.909 


178.646 


8.721 


IGR J 18027-20 16"' 


Bllb 


139.6 


4.570 


<0.2 


68 ± 1 




16.16 ±0.71 


LS 5039"- 


ON6.5 V(f) 




3.906 


0.337 ± 0.036 




236.0 ± 5.8 




XTE J 1855-026*3 


BOIaep 


361.0 


6.072 


0.04 ± 0.02 


80.5 ± 1.4 


226 ± 15 


15.19 ±0.79 


4U 1907+09'*'' 


08.5 lab 


437.5 


8.375 


0.28 ± 0.04 


83 ±2 


330 ±7 


8.75 ± 0.63 



Note. We include all HMXBs in the Galaxy and Magellanic Clouds that we have been able to find reference to in the literature for which the orbital period 
and eccentricity are reliably known. If the parameter is known to better than 3 decimal places, it is truncated and the error is smaller than 0.0 01. References 
for qu antities presented in the Table are given here f or each indiv idual system: 'Kaspi et al. 1 1994); -Harmanec et al. 12000); - Schurch et al. (2009); 

fsend e;t alj J201 1|); ^Raichur & Paul i2010h|); lAragona et al. (2009); 'Raichur & Paul (2010b); ^ Delgado-Martf et al. (2001); ''Barsukova et al.i i2005l) ; 
er. Wilson & Hagedon (1994); "Hutchines et al. ( 1985); '^Kaaret. Cusumano & Sacco (2000); ''Wilson et al. (1997); 
g, Johnston & Manchester (2004); '"Coe et al. (2007); ' ^Staubert et al. (201 1); ' '.Raichiir & PauU2010b.); "^.in't Zand. Corbet & Marsh alil JiooH); 

'-Kellev. RapDanort & Avash ( 1983); -"^Bavkal et al. (20100; ^ IScott et al.l <1997ll.lShaw et alj j2009l); ^^Finger et alj(ll999h; ^teaU owav. Wang & Morgat] 
|2005i); ^"jtVilson et al. (2003); ^' Gallowav. Morgan & LevS hm4); ^''Wilson, Finger & Camero-Arrand J2008(); ^iBaykal et all ('2007); ^ ^.Raichur & Paull 
koiOah; ^fc rundstrom et al. (200 71); 3lLevineHl99ll).lvan" der Meer et al. (2007); 3 ' Krevkenbohm et al. (2008); 3^ Raichur & Paul (2010ij 



Ray & Chakrabarty (2002); 3^ Koh et al.' ('l997^; 33ciark ( 2000); 3f Thompson et al. (2006); 3' Cusumano et al. (2010); 31^ Chakrabarty ( 1993"): 
Thompson et al. ( 2007); "' Stairs eLak (2001) ; '" Hill et al. (2005) ; "^ Aragona et al. (2009) ; "" Corbet & Mukai ( 2002) ; "^ Bavkal et al. ( 2006) . 
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the updated light curve (including over 4 years of new data) re- 
fines this period to 1 17 ± 1 d. The X-ray period seems to be driven 
by the time between outbursts, although considered with the 114 d 
optically derive d period and the ~ 620 d disk relaxation timescale 
dMcGowan et a l. 2008), a superorbital period cannot be ruled out. 
An additional argument for the 21.9 d period being orbital in na- 
ture is the re lationship between the orbital and s uperor bital peri- 
ods found bv iRaioelimanana^ Charles & Udalskil feoilh (see their 
Fig. 44). The strong correlation between these parameters suggests 
that for a superorbital period of ~ 620 d, the expected orbital pe- 
riod is around 15-20 d. During their study, those authors conclude 
that the orbital period of SXP6.85 is llOd, making the system the 
largest outlier on their plot. Our new orbital period removes this 
outlier, further improving their correlation. SXP8.80 also shows ev- 
idence for many different periodicities in it's optic al and X-ray ligh t 
curves. The possible period of 185 d reported in ICoe et al.l (|2005|) 
was only found in the MACHO red data. Those authors comment 
that they could not find the same period in the early OGLE data 
implying that the result should be taken with some caution. It is 
also very close to half a year - a period that always needs to be 
tre ated with care in such data sets. The period of 33 d reported 
by ISchmidtke & CowlevI ( |2006|) is described as 'weak' by those 
authors and possibly only seen in some early OGLE data. They 
ac cept that this could actually be the same as the 28 d reported 
bv lCorbetetai] J2004l) . Although these periodicities could be real, 
possibly a signature of stellar pulsations in the Be star, it is likely 
that they are false. Indeed, as the optical data sets have grown in size 
over the years, so has the confi dence in what is found. Hence, the 
optica l period of 28 d found bv lRaioelimanana. Charles & Udalskil 

l|) is much more likely to be real as it was found in a much 
longer data set. In this case, it is very likely to be the orbital pe- 
riod^ of SXP8.80 as it also agrees with the X-ray derived period of 
ICorbet et al. ( .2004) . 

As mentioned previously, the low metallicity environment in 
the SMC may have a significant impact on the evolution of HMXBs 
which may manifest itself as observable differences in the popula- 
tion from that of the Milky Way. Fig. |5] suggests that there is little 
difference in the relationshi p between spin and o rbital periods of 
SMC and Galactic systems. fMcBride et alj j2008l) also show there 
to be no difference in the distribution of spectral types. In the fol- 
lowing section we investigate if there is evidence of such a differ- 
ence in the distribution of eccentricities and, by inferrence, super- 
nova kick during the formation of the NS. 



4.2.1 The relationship between orbital period and eccentricity 

The orbital solutions presented here mean that for the first time, we 
can investigate the relationship between the orbital period and ec- 
centricity for a sample of SMC BeXRBs and compare this to what 
is seen in the Galaxy. In Fig. [6] we plot these two parameters for 
the 6 SMC sources and the Galactic sample. For clarity we remove 
the PSR designated systems, CI Cam and the millisecond binary 
pulsars A0538-66 and SAX J0635. 2-1-0533 from our comparison 
as these are almost certainly not 'normal' HMXBs. There are sev- 
eral things of note here. Firstly, there seems to be little difference 
between the parameter space occupied by SMC and Galactic sys- 
tems. This may suggest that the metallicity difference between the 
galaxies is not a large factor in binary evolution, although many 
more orbits of SMC systems need to be solved before we can draw 
a firm conclusion. Secondly, the Be and SG systems occupy sep- 
arate regions of the plot much like on the Corbet diagram. This is 
what we might expect considering the tight, circular orbits in SG 
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Figure 5. Corbet Diagram for all the HMXBs in Table|2]that have a known 
spin period (excluding the PSR systems). The red diamonds represent the 
SG systems, the blue triangles are the Galactic BeXRBs and the green stars 
are the SMC BeXRBs. 



systems and the wide, eccentric orbits in Be systems. The upper- 
right most red diamond (GX 301-2) has shown evide nce of be- 
havin g like a Be system despite its B I la classification jKoh et al] 
and so may be an exception to this relationship. Thirdly, 
there is a possible correlation between the two parameters, partic- 
ularly when we ex clude the small c lass of low eccentricity, long 
orbit OB transients jPfahl et al.ll2002l) in the shaded regiorj^ To try 
and quantify any possible relationship between the two parameters, 
a Spearman rank correlation coefficient was derived. This was pre- 
ferred over a linear Pearson coiTelation as the relationship seems 
to be linear-log as opposed to simply linear. The rank correlation 
coefficient for the whole dataset was computed to be 0.412 with 
a p-value of 0.007 (significant at the 3 cr level). A second corre- 
lation test was performed excluding the 8 sources proposed to be- 
long to the low eccentricity group of OB transients, resulting in a 
coefficient of 0.79 and p-value of 2.6 x 10"^ (significant at the 6(T 
level). The more significant correlation resulting from the removal 
of this small group of systems may reinforce the idea put forward 
bv lPfahl etali d2002l) that these are a separate population of Be sys- 
tems that receive a smaller supernova kick compared to 'classical' 
Be systems. It may also suggest that orbital period and eccentric- 
ity in HMXBs are somehow intrinsically related. However, fitting 
any function to the current data is difficult given the spread of val- 
ues and the unknown contribution to the trend from factors such 
as mass, orbital size and tidal circularisation timescales. Whilst the 
correlation seems significant, there is an obvious alternative to this 
scenario. Separate correlation tests on the SG systems and Be sys- 
tems yield lower correlation coefficients of 0.54 (p-value of 0.03) 
and 0.49 (p-value of 0.04) respectively. This means much of the 
overall correlation could be a result of the two groups occupying 
different parameter space. With the low eccentricity transients in- 
cluded in the Be group, this correlation coefficient drops to 0.08 (p- 



3 iPfahl et alj )2002l) argue that these objects make up a new group of 
HMXB that received a much smaller lack from the supernova explo- 
sion that created the X-ray binary. This work was originally based on 6 
systems, though now there are cu rren tly 8 members of t his group. See 
iGallowav. Wang & Morgmil )2005l) and lBavkal eTal] j2010l) for details on 
the other two systems. 
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Figure 6. Orbital period against eccentricity for all the systems presented in 
Table |2] excluding the millisecond pulsars and CI Cam. The red diamond, 
blue triangle and green star symbols represent SG, Galactic Be and SMC 
Be systems respectively, as in Fig. [3] The red and blue arrows are upper 
limits for the eccentricity of SG and Galactic Be systems respectively. The 
1 o" errors are not plotted here for claiity, but are approximately the size of 
the data point or smaller. The 8 systems within the shaded region belo ng to 
the new class of low eccentricity, long orbit HMXB tPfahl et alj2002h . 



low ec centricities and long orbital periods proposed by IPfahl et al] 
( |2002|) . or at least lies between this class and the 'normal' pop- 
ulation. On comparing the SMC sample to Galactic systems we 
find there to be little difference in their binary parameters, suggest- 
ing that the low metallicity environment in the SMC does not con- 
tribute strongly to the evolution of the binary. The different types 
of HMXB are found to occupy difi"erent regions on a plot of orbital 
period and eccentricity, similar to what is seen in the Corbet dia- 
gram. A possible correlation between the two parameters is noted, 
although it is unclear whether this is an eff'ect of the position of 
the groups on the plot or if it demonstrates a physical connection 
between the parameters in HMXB systems. The data clearly suffer 
from selection eff'ects, in particular the lack of many long orbital 
period systems that will ultimately determine the answers to the 
questions posed here. As a by-product of our work, we present a 
catalogue of the orbital parameters for every HMXB in the Galaxy 
and Magellanic Clouds for which they are known. 
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value 0.69). Again this favours the interpretation that we are seeing 
3 different populations with different evolutionary histories, rather 
than one continuous population. Conversely, we can think of the 
weak coiTelations shown in both groups as evidence for a relation- 
ship between the parameters, which is only enhanced when the two 
groups are combined. This could be telling us something signifi- 
cant about the pre-supemova phase of the binary evolutionary path. 
[Delgado-Martf et al. (2001) simulate binary periods and eccentrici- 
ties resulting from SN kicks imparted onto the NS with varying ini- 
tial orbital separations. Somewhat intuitively, the results show that 
as the initial separation is increased, a higher fraction of systems 
are found to have larger orbital periods and larger eccentricities. If 
we consider each binary system to have the same evolutionary path, 
the correlation in Fig.|6]may be telling us that SG systems gener- 
ally have tighter pre-supemova orbits than Be systems. As with the 
question of whether there is a difl"erence between SMC and galactic 
systems in this parameter space, any potential correlation in these 
data can only be confirmed by including more binary systems in the 
comparison, particularly more long orbital period Galactic systems 
and both long and short orbital period SMC systems. However, this 
will be difficult given the nature of the analysis method. Finally, we 
note the location of SXP2.37, immediately left of the shaded region 
in Fig. |6] We postulate that this source is the latest member of this 
growing class of low eccentricity OB transient HMXB and the first 
member from outside the Galaxy. 



5 SUMMARY 

The orbital parameters of four BeXRB systems in the SMC have 
been derived, taking the total number of extra-Galactic BeXRBs 
with orbital solutions to six. We find one new orbital period and 
confirm a second. Two further fits were made possible by freez- 
ing the orbital period parameter to known values from analysis of 
optical light curves. We find that one system (SXP2.37 = SMC X- 
2) is consistent with being part of the small class of HMXBs with 
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